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ABSTRACT 

Radio observations discovered large scale non thermal sources in the central Mpc 
regions of dynamically disturbed galaxy clusters (radio halos). The morphological and 
, spectral properties of these sources suggest that the emitting electrons are accelerated 

by spatially distributed and gentle mechanisms, providing some indirect evidence for 
turbulent acceleration in the inter-galactic-medium (IGM). 

Only deep upper limits to the energy associated with relativistic protons in the IGM 
have been recently obtained through gamma and radio observations. Yet these pro- 
tons should be (theoretically) the main non-thermal particle component in the IGM 
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implying the unavoidable production, at some level, of secondary particles that may 
have a deep impact on the gamma ray and radio properties of galaxy clysters. 

Following Brunetti & Lazarian (2007), in this paper we consider the advances in 
the theory of MHD turbulence to develop a comprehensive picture of turbulence in 
the IGM and extend our previous calculations of particle acceleration by compressible 
MHD turbulence by considering self-consistently the reacceleration of both primary 
and secondary particles. Under these conditions we expect that radio to gamma ray 
emission is generated from galaxy clusters with a complex spectrum that depends on 
the dynamics of the thermal gas and Dark Matter. The non-thermal emission results in 
very good agreement with radio observations and with present constraints from hard 
, X-ray and gamma ray observations. In our model giant radio halos are generated in 

mering (turbulent) clusters only. However, in case secondaries dominate the electron 
component in the IGM, we expect that the level of the Mpc-scale synchrotron emission 
in more relaxed clusters is already close to that of the radio upper limits derived by 
present observations of clusters without radio halos. Important constraints on cluster 
physics from future observations with present and future telescopes are also discussed. 
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1 INTRODUCTION 

Radio observations of galaxy clusters prove the presence 
of non-thermal components, magnetic fields and relativistic 
particles, mixed with the hot Inter-Galactic-Medium (IGM) 
(e.g. Ferrari et al, 2008). 

Potentially, cluster mergers can be responsible for the 
origin of the non-thermal components in IGM. During these 
events a fraction of the gravitational binding-energy of Dark 
Matter halos that is converted into internal energy of the 
barionic matter can be channelled into the amplification of 
the magnetic fields (e.g. Dolag et al. 2002; Subramanian et 



al. 2006; Ryu et al. 2008) and into the acceleration of par- 
ticles via shocks and turbulence (e.g. Ensslin et al 1998; 
Sarazin 1999; Blasi 2001; Brunetti et al. 2001, 2004; Pet- 
rosian 2001; Miniati et al. 2001; Fujita et al 2003; Gabici 
& Blasi 2003; Ryu et al. 2003; Cassano & Brunetti 2005; 
Hoeft & Bruggen 2007; Brunetti & Lazarian 2007; Pfrom- 
mer 2008). 

Theoretically relativistic protons are expected to be the 
dominant non-thermal particles component since they have 
long life-times and remain confined within galaxy clusters for 
a Hubble time (Volk et al. 1996; Berezinsky, Blasi & Ptuskin 
1997; Ensslin et al 1998). Proton-proton collisions in the 
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IGM inject secondary particles, including neutral pions that 
decay into gamma rays and relativistic electrons that pro- 
duce synchrotron and inverse Compton (IC) emission. So far 
only upper limits to the gamma ray emission from galaxy 
clusters have been obtained by FERMI and Cherenkov tele- 
scopes (eg., Aharonian et al 2009a, b; Aleksic et al 2010; 
Ackermann et al 2010). These upper limits, together with 
constraints from complementary approaches based on radio 
observations (eg., Reimer et al 2004; Brunetti et al. 2007) 
suggest that relativistic protons contribute to less than a 
few percent of the energy of the IGM, at least in the central 
Mpc-sized regions. 

Relativistic electrons in the IGM are nowadays studied 
by radio observations of diffuse synchrotron radiation from 
galaxy clusters. Radio halos are the most spectacular ex- 
amples of cluster-scale radio sources, they are diffuse radio 
sources that extend on Mpc-scales in the cluster central re- 
gions and are found in about 1/3 of massive galaxy clusters 
(eg. Feretti 2002; Ferrari et al. 2008; Cassano 2009). 
The origin of relativistic electrons in radio halos is still de- 
bated. In the context of the hadronic model (Dennison 1980; 
Blasi & Colafrancesco 1999; Pfrommer & Ensslin 2004) ra- 
dio halos are due to synchrotron emission from secondary 
electrons generated by p-p collisions, in which case clusters 
are (unavoidably) gamma ray emitters due to the decay of 
the 7r° produced by the same collisions. The very recent 
non-detections of nearby galaxy clusters at GeV energies by 
FERMI significantly constrain the role of secondary elec- 
trons in the non-thermal emission (Ackermann et al 2010). 
Similarly, previous works pointed out that the spectral and 
morphological properties of a number of radio halos ap- 
pear inconsistent with a simple hadronic origin of the emit- 
ting particles (eg. Brunetti et al 2008, 2009; Donnert et al 
2010a,b; Macario et al 2010). 

A second scenario proposed for the origin of radio halos is 
based on turbulent reacceleration of relativistic particles in 
connection with cluster-mergers events (eg., Brunetti et al. 
2001; Petrosian 2001; Berrington & Dermer 2003; Fujita et 
al 2003; Cassano & Brunetti 2005). The acceleration of ther- 
mal electrons to relativistic energies by MHD turbulence in 
the IGM faces serious drawbacks based on energy arguments 
(eg., Petrosian & East 2008), consequently in these models 
it must be assumed a pre-existing population of relativis- 
tic electrons in the cluster volume that provides the seed 
particles to reaccelerate during mergers. These seeds may 
be primary electrons injected by SN, AGN, galaxies and 
shocks in the cluster volume that can be accumulated for 
a few Gyrs at energies of a few hundred MeV (eg., Sarazin 
1999; Brunetti et al 2001). 

MHD turbulence theory seriously advanced in the last 
decades also affecting our view of particle acceleration in as- 
trophysical plasmas (e.g. Chandran 2000; Yan & Lazarian 
2002; Lazarian 2006a and ref therein). For this reason in 
Brunetti & Lazarian (2007) we considered the advances in 
the theory of MHD turbulence to develop a comprehensive 
picture of turbulence in the IGM and to study the reac- 
celeration of relativistic particles. In this respect, our main 
conclusions were that the compressible MHD turbulence (es- 
sentially fact modes) , generated in connection with energetic 
cluster-mergers, is the most important source of stochastic 
particle reacceleration in the IGM and that the interaction 
between this turbulence and the relativistic electrons may 



explain the origin of radio halos, provided that enough seed 
relativistic electrons are accumulated in the IGM. 
At the same time relativistic protons should be accumu- 
lated (at some level) in the IGM and consequently MHD 
turbulence may also reaccelerate these protons and their 
secondary products (Brunetti & Blasi 2005). On one hand, 
in this case secondaries come into the picture as a natural 
pool of seed particles, on the other hand the reacceleration 
of relativistic protons is expected to enhance the produc- 
tion rate of these secondaries with the possible drawback 
that this would overcome the effect of reacceleration, thus 
making the overall picture unsustainable. 
Previous papers addressed this problem and model turbulent 
reacceleration of both primary and secondary particles in 
the IGM in the context of cluster-mergers driven turbulence 
(Brunetti & Blasi 2005; Brunetti et al 2009a). These papers 
explored the complex non-thermal spectrum from galaxy 
clusters that is expected in this situation and first predicted 
the interplay of a long-living spectral component, due to 
relativistic protons and their chain of secondary particles, 
with a transient emission, leading to the formation of radio 
halos, that is due to the reacceleration of electrons in clus- 
ter mergers. Remarkably, in this general situation gamma 
ray emission is expected from galaxy clusters, although at a 
lower level than that expected from the pure hadronic case 
(eg. Brunetti et al 2009a)*. 

All these papers addressing the reacceleration of primary 
and secondary particles focus on Alfvenic reacceleration and 
assume the injection of Alfven modes at small, resonant, 
scales, in which case however it is difficult to derive a overall 
self-consistent picture connecting clusters mergers and the 
generation of these modes at such small scales. This makes 
necessary an exploration of this picture in the context of a 
more comprehensive model of MHD turbulence in clusters, 
that is the goal of our paper. 

Our study is also timely due to the recent gamma and radio 
observations that put severe constraints to the energy den- 
sity of the relativistic proton component in galaxy clusters 
(Brunetti et al 2007; Aharonian et al 2009a,b; Aleksic et al 
2010; Ackermann et al 2010), and that allows for including 
this component in models with substantially less degree of 
freedom than in the past. 

For these reasons, in this paper we extend our calculations 
presented in Brunetti & Lazarian (2007) by considering self- 
consistently relativistic protons, the generation of secondary 
particles in the IGM and their reacceleration by MHD tur- 
bulence generated in cluster mergers. 

In Section 2 we discuss the properties of turbulence in 
galaxy clusters, in Sect. 3 we summarize the formalism for 
particle acceleration and evolution, in Sect. 4 we first discuss 
the connection between mergers and particle reacceleration 
by turbulence and then show our results, and in Sect. 5 we 
give our conclusions and discuss model simplifications, fu- 
ture extensions of the work and a comparison with previous 
works. 

A ACDM cosmology with H = 70 km/s/Mpc, A = 0.7 and 
Ho = 0.3 is assumed. 



* assuming the same cluster magnetic field and for a given syn- 
chrotron luminosity 
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2 TURBULENCE IN THE IGM 

Numerical simulations of galaxy clusters suggest that tur- 
bulent motions may store an appreciable fraction, 5-30%, of 
the thermal energy of the IGM (e.g., Roettiger, Burns, Lo- 
ken 1996; Roettiger, Loken, Burns 1997; Ricker & Sarazin 
2001; Sunyaev, Bryan & Norman 2003; Dolag et al. 2005; 
Vazza et al. 2006, 2009a; Paul et al 2010). The largest tur- 
bulent eddies should decay into a turbulent velocity field on 
smaller scales, possibly developing a turbulent cascade. 

It is known that the mean free path of thermal pro- 
tons arising from Coulomb collisions in the hot IGM may 
be very large, ten to hundred kpc. Fluids in such a colli- 
sionless regime can be very different from their collisional 
counterparts (Schekochihin et al. 2005; 2010). The parallel 
to magnetic field viscosity of IGM can be very large and 
the collisionless plasmas are subject to various instabilities. 
Those, however, we believe change the effective collisionality 
of the fluid, justifying the application of MHD to describing 
the IGM at least on its large scales (see Lazarian et al. 2010). 
Indeed, particles in plasmas can interact through the media- 
tion of the perturbed magnetic fields, and effective collision- 
ality of plasmas may differ dramatically from the textbook 
estimates. The difference stems from the very instabilities 
that are present in the IGM plasmas (firehose, mirror, gy- 
roresonance etc.). These instabilities are expected to transfer 
the energy from the turbulent compressions on the scales less 
or equal to the particle mean free path to the perturbations 
at the particle gyroscale. This has been described in Lazar- 
ian & Beresnyak (2006), as a result of the scattering, the 
mean free path of particles decreases, which makes the fluid 
essentially collisional over a wide range of scales, with the 
critical scale for which the fluid gets effectively collisional 
that is expected to decrease with the increase of turbulent 
driving rate. While some parts of the aforementioned pa- 
per dealing with the interaction of turbulence and cosmic 
rays remain controversial, a similar approach should be re- 
liably applicable to thermal plasma particles. Therefore we 
believe that on scales much larger that the thermal parti- 
cle gyroradius the turbulence can be treated in the MHD 
approximation and shall apply below the theory of MHD 
turbulence to the IGM. 

Turbulence generated during cluster mergers is ex- 
pected to be injected at large scales, L ~ 100 — 400 kpc, 
with typical velocity of the turbulent eddies at the injec- 
tion scale around V ~ 300 — 700 km/s (eg. Subrama- 
nian et al 2006). This makes turbulence sub-sonic, with 
M a = Vo/c a « 0.25 — 0.6, but strongly super- Alfvenic, 
with Ma = Vo/va ~ 5 — 10. Turbulent motions at large 
scales are thus essentially hydrodynamics and the cascading 
of compressive (magnetosonic) modes may couple with that 
of solenoidal motions (Kolmogorov eddies) . 

In Brunetti & Lazarian (2007) we discussed that the im- 
portant consequence of the turbulence in magnetized plasma 
is that both solenoidal and compressive modes in hot galaxy 
clusters would not be strongly affected by viscosity at large 
scales and an inertial range is established, provided that the 
velocity of the eddies at large scales exceeds « 300 km/s. 
In the Kolmogorov cascade the turbulent velocity Vi scales 
as Vl(Z/L ) 1//3 , and at scales less than Ia ~ L M^ 3 the 
turbulence gets sub-Alfvenic and we enter into the MHD 
regime (see discussion in Lazarian 2006b). For the parame- 



ters above the scale I a ~ 0.1 — 1 kpc, but the actual number 
of Ma is not certain and therefore our estimate of Ia should 
be treated with caution. Fortunately, the above uncertainty 
does not change the results of our paper appreciably. 

Compressible MHD turbulence is a subject where a 
number of important insights have been obtained much be- 
fore these ideas can be tested; the pioneering works in the 
area include Montgomery & Turner (1981), Shebalin et al. 
(1983), Higdon (1984). The key idea of critical balance by 
Goldreich & Sridhar (1995) has influenced in a profound way 
our further thinking of the MHD cascade. 
In the MHD regime, at smaller scales where Vi < va, three 
types of modes should exist in a compressible magnetized 
plasma: Alfven, slow and fast modes. Slow and fast modes 
may be roughly thought as the MHD counterpart of the 
compressible modes, while Alfven modes may be thought as 
the MHD counterpart of solenoidal Kolmogorov eddies (a 
more extended discussion can be found in Cho, Lazarian & 
Vishniac 2002, and ref. therein). 

Turbulence in the IGM is most likely a complex mixture 
of several turbulent modes. We shall assume that a sizeable 
part of turbulence at large scales (namely at scales where 
the magnetic tension does not affect the turbulent motions) 
is in the form of compressible motions. This is reasonable as 
these modes are easily generated in high beta medium (eg. 
Brunetti & Lazarian 2007 and ref. therein). 
Situation may be radically different at smaller scales where 
the magnetic field tension affects turbulent motions, i.e. in 
the MHD- regime, I < I a - In this case, MHD numerical sim- 
ulations have shown that a solenoidal turbulent forcing gets 
the ratio between the amplitude of Alfven and fast modes 
in the form (Cho & Lazarian 2003) : 

which essentially means that coupling between these two 
modes may be important only at I w Ia (in the MHD- 
regime it should be (8V) S < va) since the drain of energy 
from Alfvenic cascade is marginal when the amplitudes of 
perturbations become weaker. Most importantly in galaxy 
clusters it is » v\ and thus the ratio between the am- 
plitude of Alfven and fast modes at scales / < Ia is ex- 
pected to be small, {8V)l/{5V)l < {v A /c s ) 2 ~ 10~ 2 (this 
for solenoidal forcing at I « I a)- 

A more recent work by Kowal & Lazarian (2010) decom- 
posed turbulent motions into slow, fast and Alfven modes 
using wavelets. This approach is better justified than the 
decomposition in Fourier space employed in Cho & Lazar- 
ian (2002, 2003). Indeed, Alfven and slow modes are defined 
in the local system of coordinates (see Lazarian & Vishniac 
1999, Cho & Vishniac 2000, Maron & Goldreich 2001) and 
therefore the procedure of Fourier decomposition by Cho & 
Lazarian (2002, 2003, see 2005 for a review) can only be sta- 
tistically true. The wavelet decomposition procedure is more 
localized in space and therefore is potentially more precise. 
Nevertheless, the results on mode decomposition in Kowal 
& Lazarian (2010) agree well with those in Cho & Lazarian 
(2003), which provides us with more confidence about the 
properties of MHD turbulence that we employ in the paper 
to describe turbulence interaction with energetic particles. 
Our treatment of turbulence assumes that the turbu- 
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lence is balanced, i.e. the energy flux of waves^ moving in 
one direction is equal to the energy flux in the opposite di- 
rection. While the theory of imbalanced turbulence is being 
intensively developed (see Lithwick & Goldreich 2001; Beres- 
nyak & Lazarian 2008, 2009, 2010; Chandran 2008, Perez & 
Boldyrev 2009) we do not expect that the effects of imbal- 
ance would dominate in the cluster environments. First of 
all, the turbulence driving is not expected to be strongly 
localized and then the effects of compressibility should de- 
crease the local turbulence imbalance. 



dp 2 



N P (p,t) 
r pp (p) 



(3) 



where \dp/dU\ marks Coulomb losses (Sect. 3.3), D pp is the 
diffusion coefficient in the momentum space of protons due 
to the coupling with magnetosonic modes (Sect. 3.2), and 
t pp is the proton life-time due to pp collisions in the IGM 
(Sect. 3.3; see also Ensslin et al 2007). 

For isotropic turbulence (Sect. 3.2) the diffusion equation in 
the k-space is given by : 



3 PARTICLE ACCELERATION, ENERGY 
LOSSES AND SECONDARY PARTICLES 

In this paper we provide an extension of our previous cal- 
culations. We assume the picture of MHD turbulence in the 
IGM, as derived in Brunetti & Lazarian (2007), to calcu- 
late the reacceleration of relativistic particles by compress- 
ible MHD turbulence by taking into account self-consistently 
also the generation and reacceleration of secondary particles. 
The aim of this Section is to present the formalism and the 
main assumptions used in our calculations. 



3.1 Basic formalism 

We model the re-acceleration of relativistic particles by 
MHD turbulence in the most simple situation in which only 
relativistic protons are initially present in a turbulent IGM. 
These protons generate secondary electrons via p-p collisions 
and in turns secondaries (as well as protons) are reacceler- 
ated by MHD turbulence. 

We model the time evolution of the spectral energy 
distribution of electrons, A^ e _ , and positrons, AT+, with an 
isotropic Fokker-Planck equation^ : 



dt 



d_ 

dp 



N; 



dp 

~dl, 



dp 


)} 


~dli 





dp 2 

+Qf[p,t;N p (p,t)], 



[D±N±(p,t)] 



(2) 



where \dp/dt\ marks radiative (r) and Coulomb (i) 
(Sect. 3.3), D pp is the electron/positron diffusion coefficient 
in the momentum space due to the coupling with magne- 
tosonic modes (Sect. 3.2), and the term Q^r accounts for the 
injection rate of secondary electrons and positrons due to 
p-p collisions in the IGM (Sect. 3.4). 

The time evolution of the spectral energy distribution of 
protons, N p , is given by : 



dN p (p,t) = 8_ 
dt dp 



t Strong MHD turbulence presents a case of dualism of waves and 
eddies (see Lazarian & Vishniac 1999; Cho, Lazarian & Vishniac 
2003). Strong non-linear damping of oppositely moving Alfvcnic 
wave packets makes them act like eddies, 
i isotropization of particles is discussed in Sect. 3.2. 



dt 



■5^ri(M)W(M) , 



(4) 



where D kk is the diffusion coefficient in the k-space, Ti(k, i) 
are the different damping terms, and I(k, t) — I S(k — k ) is 
the turbulence injection term, i.e. we consider the most sim- 
ple situation where turbulence is injected at a single scale, 
with wavenumber k . The wave-wave diffusion coefficient 
of magnetosonic modes (Kraichnan treatment) is given by 
(Brunetti & Lazarian 2007 and ref therein): 



D k 



{V ph )k' 



( W(M) 



(5) 



where (V p h) is a representative, averaged (with respect to 
8), phase velocity. 

We shall assume isotropic MHD turbulence to calcu- 
late the particle acceleration rate at any time. This assump- 
tion is appropriate for super-Alfvenic turbulence and fast 
modes (e.g., Cho & Lazarian 2003), provided that collision- 
less dampings are not efficient. At smaller scales, collisionless 
dampings with thermal particles in the IGM become severe 
and modify the spectrum of turbulent modes, and a cut-off 
in the turbulent spectrum is generated at k — k c , where 
the damping time-scale becomes shorter than the cascad- 
ing time. The most important damping with thermal (and 
relativistic) particles is the Transit-Time-Damping (TTD) 
that is highly anisotropic (e.g. Schlickeiser & Miller 1998; 
Brunetti & Lazarian 2007; Yan et al. 2008) being stronger for 
6 ~ 7r/2 and causing the spectrum of the turbulent modes to 
become anisotropic at scales k ~ k c (0). On the other hand, 
in Brunetti & Lazarian (2007) we have shown that, under 
physical conditions typical of the IGM, hydro-motions bend 
the magnetic-field lines in a time-scale, Tbb ~ Ia/va, that 
is comparable to the damping time-scale of turbulence, in 
which case there is chance that approximate isotropization 
of the turbulent spectrum is maintained even at scales where 
dampings are severe. 

Thus following Brunetti & Lazarian (2007) we shall assume 
a simplified turbulent (isotropic) spectrum in the form : 



W(fc)« (I oP (V pb ))i k~ 



(6) 



for k < k < fee, with the cut-off wavenumber estimated 
from the condition that the damping time-scale becomes 
smaller than the cascading time scale : 



Td « 1/(^2 Ti{k, 6)} = r kk » k 2 /D k 
that is : 



(7) 
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h c — C 



(8) 



where quantities (..} are averaged with respect to 6 and ^ I\ 
is the total (thermal and non-thermal) damping term due 
to TTD resonance and the constant C ~ a few (Brunetti & 
Lazarian 2007; see also Matthaeus & Zhou 1989, for details 
on Kraichnan constants). 



3.2 Momentum diffusion coefficient due to 
compressible MHD turbulence 

Compressible turbulence can affect particle motion through 
the action of the mode-electric field via gyroresonant inter- 
action (e.g., Melrose 1968), the condition for which is : 



k \\ v \\- n -=° 



(9) 



where n = ±1, ±2, .. gives the first (fundamental), second, 
.. harmonics of the resonance, while v» = /j,v and k\\ = r\k 
are the parallel (projected along the magnetic field) speed 
of the particles and the wave-number, respectively. 

Following Brunetti & Lazarian (2007) we assume 
particle-mode coupling through the Transit-Time Damp- 
ing (TTD), n = 0, resonance (e.g., Fisk 1976; Eilek 1979; 
Miller, Larosa & Moore 1996; Schlickeiser & Miller 1998). 
In principle, this resonance changes only the component of 
the particle momentum parallel to the seed magnetic field 
and this would cause an increasing degree of anisotropy 
of the particle distribution leading to a less and less effi- 
cient process with time. Thus an important aspect in this 
working picture is the need of isotropization of particle mo- 
menta during acceleration (e.g., Schlickeiser & Miller 1998). 
In this paper we shall assume continuous isotropization of 
particle momenta. Isotropy may be provided by several pro- 
cesses discussed in the literature. These include electron 
firehose instability (Pilipp & Volk 1971; Paesold & Benz 
1999), and gyro-resonance by Alfven (and slow) modes at 
small scales, provided that these modes are not too much 
anisotropic (cf. Yan & Lazarian 2004) § . Gyro-resonance may 
also occur with electrostatic lower hybrid modes generated 
by anomalous Doppler resonance instability due to pitch an- 
gle anisotropies (e.g., Liu & Mok 1977; Moghaddam-Taaheri 
et al. 1985) and, possibly, with whistlers (e.g., Steinacker 
& Miller 1992). In addition, Lazarian & Beresnyak (2006) 
proposed isotropization of cosmic rays due to gyroresonance 
instability that arises as the distribution of cosmic rays gets 
anisotropic in phase space. 

We adopt the momentum-diffusion coefficient of parti- 
cles, D pp , as derived from detailed balancing argument, i.e. 
relating the diffusion coefficient of a a-species to the damp- 
ing rate of the modes themselves with the same particles 
(eg., Eilek 1979; Brunetti & Lazarian 2007)^ : 



n r\ %2 2 1 r \ m r2 sin 3 (fl) „„., V ph /c. 



V ph /c 



cos(0)| 

2 



cos 9 
dkW B (k)k (10) 



where B is the background (unperturbed) magnetic field, 
k c is given in Eq. 8 and 



^B(k) = ±( ^ l [ B J* \W(k) , 



(11) 



I3 pl " WnW(k) 

where \Bk\ 2 /VV(fc) is the ratio between magnetic field fluc- 
tuations and total energy in the mode (the quantity (..) in 
Eq. 11 indicates average with respect to 9 and is of order 
unity, see Brunetti & Lazarian 2007 for further details). 

3.3 Energy losses for electrons and protons 

The energy losses of relativistic electrons in the IGM are 
dominated by ionization and Coulomb losses, at low ener- 
gies, and by synchrotron and IC losses, at higher energies 
(eg. Sarazin 1999). The rate of losses due to the combina- 
tion of ionization and Coulomb scattering is (in cgs units): 



(dp\ 
\dtJi 



-3.3 x 10"^n th 



1 + 



ln(7/w t h) 
75 



(12) 



where nth is the number density of the thermal plasma. The 
rate of synchrotron and IC losses is (in cgs units): 

2 „;„2 a 



(±) =_ 

V dt ) rad 



4.8 x 10 



-4 2 

V 



my 



2/3 



(13) 



where B p g is the magnetic field strength in units of /j,G, 
and 9 is the pitch angle of the emitting leptons; in case of 
efficient isotropization of the electron momenta, the sin 2 9 is 
averaged to 2/3. 

For relativistic protons, the main channel of energy 
losses in the IGM is provided by inelastic p-p collisions. The 
life-time of protons due to pp collisions is given by : 

r pp (p) = 1 (14) 

In this paper we use the inclusive cross section, tj ± '°(p p ), 
given by the fitting formulae in Dermer (1986a) which allow 
to describe separately the rates of generation of neutral and 
charged pions. 

For trans-relativistic and mildly relativistic protons, energy 
losses are dominated by ionization and Coulomb scattering. 
Protons more energetics than the thermal electrons, namely 
with p p > p c ee (3/2m e /m p ) 1 / 2 /3 e (ft ~ 0.18(T/10 8 iO 1/2 
is the velocity of the thermal electrons), are affected by 

Coulomb interactions. Defining x m = (^j^j fie, one has 
(Schlickeiser, 2002): 



''^ - -1.7 xlO- 26 



\dtJi 



\ 10- 3 J 



(cgs) 



(15) 



§ The latter condition means that the Alfvenic modes are consid- 
ered for scales not much less than I a , provided that the turbulence 
injection is isotropic 

K This is obtained by assuming that particles isotropy is main- 
tained during acceleration (see Brunetti & Lazarian 2007) 



3.4 Injection of Secondary Electrons 

The decay chain that we consider for the injection of sec- 
ondary particles in the IGM due to p-p collisions is (Blasi 
& Colafrancesco 1999): 
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p + p — > 7T° + 7r + + 7r + anything 



(16), (17) and (21), we obtain the rate of production of sec- 
ondary electrons/positrons : 



it — > 77 



that is a threshold reaction that requires protons with ki- 
netic energy larger than T p « 300 MeV. 

A practical and useful approach to describe the pion 
spectrum both in the high energy (E p > 10 GcV) and 
low energy regimes was proposed in Dermer (1986b) and 
reviewed by Moskalenko & Strong (1998) and Brunetti & 
Blasi (2005), and is based on the combination of the iso- 
baric model (Stecker 1970) and scaling model (Badhwar et 
al., 1977; Stephens & Badhwar 1981). 
The injection rate of pions is given by : 

Q±-°(E,t) = n* h c ( dpN p (p,t)p /^ Ep)a± '°^\ (16) 
Jp, V 1+ \ m P c IPp) 

where we adopted F n as given in Brunetti & Blasi (2005) 
that use the isobaric model for E p < 3 GeV, a scaling model 
for E p > 10 GeV and a linear combination of the two models 
for intermediate energies, and where ptr is the threshold 
momentum of protons for the process to occur. 

The injection rate of relativistic electrons/positrons is 
given by : 



Qe±(p,t) = J Qn(E lr ±,t)dE 7r j dE^ x 

F e ±(E„,E^,E e )F^E^E„), (17) 

where F^(E e , E^, En) is the spectrum of electrons and 
positrons from the decay of a muon of energy E^ produced 
in the decay of a pion with energy E n , and F P ,(E P ,, E n ) is the 
muon spectrum generated by the decay of a pion of energy 
Ett that is : 

™2 -, 



F^(E^, En) 



(18) 



ml - ml y/E'i - ml 

between a kinematic minimum and maximum muon energy 
given by : 



(19) 



where 7^ is the Lorentz factor of the muon in the pion 
frame, fa ~ 0.2714 (Moskalenko & Strong 1998), and from 
kinematics : 



w- m 7 m Pm = 



m'j - ml 
2m 7r 



(20) 



In order to simplify calculations, following Brunetti & Blasi 
(2005), we assume that the spectrum of muons is a delta- 
function : 



Fp(Ep, E n ) = 8 [Ey, — E^ 
where 



E u = 



ml - ml Ek 
ml 2/31 



(21) 



(22) 



We use the spectrum of electrons and positrons from the 
muon decay, F e ±(E 7T , _E M , E e ), as given by Blasi & Co- 
lafrancesco (1999), and combining their results with Eqs. 



Qe±(p,t) = 



ml -ml Je 



E, 



Ett/3^ 



3 P N(E P ) x 



1- {^) 2 o ± {E p )F^,E p )F, (23) 



where E min = 2E e fam 2 r /(ml - m 2 ,), and 



F -12-4 X +3 A "2^6 



(^(^ + i)A 2 + (^ + i)A 3 ), 



for ^— &i < A < 1; 
1 + fa " " 



and 



X 2 fa 



(i - fa) 2 



3 \l-fa 



2A 2 



1- 



+ A(1 + , 



2X 3 {3 + /3 



for < A < 

and where we use the following definitions : 
T E e E*( 2/3> 2 



3(1 -fa) 2 

1-fa 



1 + 



and 



Pu = ~ 



fa (ml - ml) 2 



4^ -l + (^) 4 



(24) 
(25) 

(26) 



4 RESULTS 

In this Section we calculate particle acceleration and non- 
thermal emission in galaxy clusters by assuming that MHD 
turbulence is generated during mergers between clusters and 
reaccelerates primary and secondary particles. 
Following previous works (eg. Brunetti & Blasi 2005) we 
adopt a simplified situation where we do not consider pri- 
mary electrons in the IGM (see Sect. 5 for discussion), and 
where the two main ingredients are (i) relativistic protons, 
that are believed to be the most important non-thermal 
particle components in the IGM (eg., Blasi et al 2007 for 
review), and (ii) the MHD turbulence. On one hand, rela- 
tivistic protons inject secondary particles via p-p collisions 
in the IGM that produce radiation from the radio to the 
gamma ray band, at the same time MHD turbulence may 
reaccelerate relativistic protons and secondary electrons in 
the IGM, generating radio halos and leaving an imprint in 
the general non-thermal properties of galaxy clusters. 

At variance with the aforementioned paper that focus 
on the Alfvenic case, following the model of MHD turbu- 
lence in galaxy clusters by Brunetti & Lazarian (2007), wc 
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Figure 1. Left panel: The life— time of relativistic electrons in the IGM at z = 0.2 as a function of their Lorentz factor. Thick (blue) 
lines are for cluster cores (B = 3/iG, n t ^ = 2 X 10 — 3 cm -3 ) and thin (red) lines are for cluster periphery (B = 0.5^tG, n t ^ = 10 — 4 cm -3 ). 
We report the total life-time (solid lines) and the life-times due to single processes: Coulomb losses (dashed lines), synchrotron and 
IC losses (dotted lines), and bremsstrahlung losses (long dashed lines). The yellow region marks the turbulence eddy turnover time 
(assuming L ~ 200 - 300 kpc, (V L /c s ) 2 ta 0.1 - 0.3, T « 10 s K), while the blue region marks the range of Lorentz factors of the 
relativistic electrons emitting at frequencies RJ 300 — 1400 MHz. Right panel: The life-time of cosmic— ray protons in the IGM as a 
function of the particle momentum. Thick (blue) lines are for cluster cores and thin (red) lines are for cluster periphery. We report the 
total life— time (solid lines) and the life-times due to single processes: Coulomb losses (dotted lines) and pp-collisions (dashed lines). The 
yellow region marks the turbulence eddy turnover time, while the blue region marks the range of momentum of relativistic protons that 
mostly contribute to the injection of secondary electrons emitting at frequencies £3 300 — 1400 MHz. 



assume that MHD turbulence in the IGM is in the form 
of compressible modes whose cascading from large to small 
scales results in a isotropic turbulent-spectrum. This allows 
us to readily connect the injection of turbulence at large 
scales with the particle acceleration process and to study 
the theoretical framework of the connection between cluster 
mergers and turbulent reacceleration of relativistic particles 
in the IGM. As a matter of fact calculations reported in this 
Section provide an extension of those in Brunetti & Lazarian 
(2007) that consider reacceleration of (only) primary parti- 
cles by compressible MHD modes. 

As already stressed, one of the main motivations for these 
new calculations comes from the recent gamma ray and ra- 
dio observations that put severe constraints on the energy 
density of relativistic protons in galaxy clusters (Brunetti et 
al 2007; Aharonian et al 2009a,b; Aleksic et al 2010; Acker- 
mann et al 2010), allowing for including secondary particles 
in turbulent -acceleration models with substantially less de- 
gree of freedom than in the past. 

4.1 Turbulent reacceleration, time— scales and 
connection with mergers 

In the framework adopted in our paper the idea is that ra- 
dio halos are generated by the reacceleration of relativis- 
tic electrons (secondaries in our specific case) by turbulence 
generated during cluster mergers. 

The leading processes in the context of our scenario are the 
generation, cascading and dissipation of turbulence in the 
IGM, and the acceleration and cooling of relativistic parti- 
cles. These processes, and the cluster-cluster collisions them- 
selves, have their own time-scale and the general picture 
breaks out from the interplay of all these time-scales. Fig- 



ure 1 reports the time-scales of the most relevant processes 
as calculated in different, relevant, regions of the cluster vol- 
ume. 

The generation and dissipation of turbulent motions in 
the IGM is not studied in great details, however present nu- 
merical simulations suggest that these motions can be gener- 
ated in galaxy clusters for a substantial fraction of the period 
of cluster-cluster interaction, a few Gyr, possibly driven by 
shock waves that cross the cluster volume and by the slosh- 
ing (and stripping) of cluster cores (e.g., Dolag et al. 2005; 
Vazza et al. 2009a; Paul et al. 2010; ZuHone et al. 2010). 
Under our working picture, the compressible turbulence, af- 
ter being injected at larger scales, decays at smaller scales 
where it dissipates through dampings with thermal and non- 
thermal particles in the IGM (Sect. 3.1). The turbulence 
decay requires about one eddy turnover time, that in the 
case of fast mode is r^k ~ ^j-(Lfk) 1 ^ 2 (e.g., Yan & Lazar- 

L 

ian 2004), implying a unavoidable delay between the first 
generation of large-scale turbulence in a given region and 
the beginning of the phase of particle reacceleration (that 
is mainly due to the non linear interaction of particles with 
turbulent modes at smaller scales) in the same region. We 
believe however that this delay does not break the temporal 
connection between mergers and particle acceleration, since 
the eddy turnover time of compressible turbulence in mas- 
sive (hot) clusters is Tkk ~ 0.2 — 1 Gyr (by assuming typical 
injection scales L a « 200-300 kpc and (V L /c s ) 2 w 0.1-0.3), 
that is smaller than the typical duration of cluster-cluster 
interaction. 

Under our simplified working picture, where turbulence is 
(only) generated by energetic cluster mergers, compress- 
ible turbulence dissipates completely in a few eddy turnover 
times, as soon as galaxy clusters becomes more relaxed. The 
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Figure 2. A scheme of the processes taken into account in our 
calculations and the coupling between them. The starting points 
are the properties of the magnetized IGM (n t j,, T, B ), the initial 
spectrum of relativistic protons (N p (p, 0)) and the injection rate 
of turbulence (I(k)). 

decrease of the efficiency of turbulent-particle acceleration 
and the suppression of non-thermal cluster-scale emission, 
i.e. the "dissipation" of radio halos, are even faster due to 
the fact that (i) the acceleration efficiency scales non lin- 
early with the turbulent spectrum 1 1 and (ii) the cooling time 
of the radio-emitting electrons is short, ~ 0.1 Gyr (Fig. 1). 
The consequence is a tight connection between radio halos 
and cluster mergers, although the picture may be more com- 
plex as the secondary particles, continuously injected in the 
IGM, should generate synchrotron emission at some level 
also in relaxed clusters (see Sect. 4.3.2). 
Radio observations of statistical samples of galaxy clusters 
show that (i) giant radio halos form (only) in merging clus- 
ters, (ii) that their life-time in merging clusters is of the 
order of 1 Gyr, and (iii) suggest that halo emission should 
dissipate in relaxed clusters in a short, < Gyr, time-scale 
(e.g. Hwang 2004; Brunetti et al. 2007, 2009b; Venturi et 
al 2008); these observational points are consistent with the 
working picture described in this Section. 

4.2 Spectral evolution of reaccelerated protons 
and electrons 

In this Section we report on some relevant results on the 
evolution of the particles spectrum (protons and secondary 
electrons/positrons) subject to TTD resonance with com- 
pressible MHD turbulence; the approach followed in this sec- 
tion is than used in next Section to calculate non-thermal 
emission from galaxy clusters. 

We adopt a simplified situation: Figure 2 shows the 
chain of physical processes that we consider and their inter- 
play. We assume that the thermal IGM is magnetized and, 
at time = 0, consider (only) relativistic protons, with initial 
spectrum N p (p) = K p p~ 26 . The presence of relativistic and 
thermal protons determines the initial efficiency of injection 

II e.g., T acc <x D pp <x W?, combining Eqs. 6, 8 and 10, see also 
Brunetti & Lazarian (2007) 



of secondary particles in the IGM (Sect. 3.4), whose initial 
spectrum is calculated assuming stationary conditions (e.g. 
Dolag & Ensslin 2000; i.e. by taking D pp = and dN/dt = 
in Eq. 2). 

Following Brunetti & Lazarian (2007), we assume that com- 
pressible turbulence is injected at large scales and develops 
a quasi-stationary spectrum (Eq. 6) due to the interplay 
between non-linear wave-wave interaction and collisionless 
dampings at smaller scales (Sect. 3.1). Following Brunetti & 
Lazarian (2007) we consider a time-independent damping, 
J^T; ~ Tth, that is obtained by assuming the (initial) phys- 
ical properties of the IGM. This is motivated by the fact that 
in our model TTD dampings with relativistic particles are 
sub-dominant with respect to those with thermal IGM (see 
Cassano & Brunetti 2005 and Brunetti & Lazarian 2007) 
and by the fact that the thermal properties of the IGM are 
not greatly modified by turbulence. 

In our calculations we follow self-consistently the re- 
acceleration of relativistic protons due to TTD with the 
compressible turbulent-modes, the generation (and its evo- 
lution with time due to the evolution of the spectrum of pro- 
tons) of secondary electrons and positrons through collisions 
between these protons and the IGM, and the reacceleration 
of the secondaries (see Fig. 2) . 

Theoretically we expect that magnetic field can be amplified 
by turbulence in the IGM (e.g. Subramanian et al 2006), al- 
though as a necessary simplification we do not include that 
amplification process in our calculations. The main reason 
is that the spectrum of the reaccelerated relativistic elec- 
trons is expected to evolve more rapidly than the magnetic 
field in the IGM (e.g. Cassano 2010). We also stress that 
present data do not show a clear connection between the 
magnetic field properties and cluster dynamics (e.g., Clarke 
et al 2001; Govoni 2006) that leaves the process of magnetic 
field amplification still poorly constrained. 

Turbulent acceleration can be thought as the combina- 
tion of a systematic effect, that causes the boosting of the 
spectrum of particles at higher energies, and a stochastic ef- 
fect, that causes a broadening of the spectrum with no net 
acceleration (e.g. Melrose 1968; Petrosian 2001). The time- 
scale of the systematic acceleration is : 

_ P 3 _ jf_ (97) 
acc dp 2 D pp /dp 4D PP V ' 

that does not depend on the particle energy in the case of 
TTD acceleration (Sect. 3.2), provided that the spectrum of 
compressible turbulence is isotropic (Sect. 3.1). 
By considering a reference value of the acceleration time due 
to TTD resonance in the IGM, r acc w 10 s yrs (e.g., Cassano 
& Brunetti 2005; Brunetti & Lazarian 2007), Coulomb** 
and radiative losses in the IGM prevent the reacceleration 
of electrons with energies E < 10 MeV and E > 10 GeV, 
respectively (Fig.l). 

On the other hand, TTD resonance in the IGM may 
reaccelerate supra-thermal protons up to high energies, and 
consequently in our model the energy density of relativistic 
protons increases with mergers, similarly to the energy of 
the thermal IGM. By considering a toy scenario where these 



** In the external regions of galaxy clusters Coulomb losses are 
less severe and lower energy electrons can be reaccelerated (Fig.l) 
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Figure 3. The time evolution of the spectrum of electrons (left panel) and protons (right panel). From bottom to top, calculations are 
shown after 0.3, 1, 2, 4, 6.5 xlO 8 yrs of rcacceleration. Calculations are obtained at z = 0, assuming (Vl/c s ) 2 = 0.22, with injection 
scale L = 300 kpc, T = 10 s K, thermal density n th = 2.64 X lO^cm" 3 and B = 4.69AtG. 



protons are simply injected in the cluster volume, with ini- 
tial energy density £.rc(0), and then reaccelerated by MHD 
turbulence during major cluster mergers, the rate of increase 
of their energy is : 



decn/dt~ I dkW{k,t)T C R{k,t) 



(28) 



where the TTD damping of compressible MHD turbulence 
by relativistic protons in case p p i >> 1 (from Eq. 31 in 
Brunetti & Lazarian 2007) is : 



Tcn/oJr 



7T 2 c a ecu sin 2 9 
4 c t t h | cost? 



d£p 

dp f 



(29) 



where {...} = s + 2 assuming a power law energy distribution 
of relativistic protons N(p) oc p~ 3 . By assuming a "sonic" 
turbulent-forcing, I t c i ~ eth, where t c i ~ 3 — 6 Gyrs is the 
life-time of massive clusters, and a typical merging history 
of massive clusters (eg., Cassano & Brunetti 2005 and ref. 
therein), we expect ecu < few percent of the thermal energy 
density, provided that relativistic protons are injected in the 
IGM with an "initial" energy density ecfl(0)/eth ~ 0.001 — 

o.oitt. 

In Figure 3 we report the time evolution of relativistic 
electrons and protons in a hot, T ~ 10 s K, IGM by assuming 
(Vl/cs) 2 =0.22, in which case the TTD acceleration time is 
Tacc ~ 10 8 yrs. Radiative losses prevent the acceleration of 
relativistic electrons above a maximum energy, "/max ~ 10 4 , 
producing a bump in the spectrum of electrons (Fig. 3a). 
The number density of high energy electrons increases with 



this is obtained by using the approximate scalings ^ - ' ss I 
in Eq.28 



time also for 7 > 'ymax (Fig- 3a), that is because the in- 
jection rate of secondary electrons is enhanced with time as 
protons are accelerated (Fig. 3b). Overall the mechanism is 
very efficient : on one hand the acceleration of relativistic 
protons enhances the injection rate of secondary electrons, 
at the same time an increasing number of secondary elec- 
trons accumulates at energies w "/max where cooling is bal- 
anced by acceleration. The combination of these two effects 
boosts the spectrum of electrons at energies 7 w "/max- 

In Figure 4 we report the evolution of the spectrum of 
electrons with time assuming physical conditions that span 
the volume occupied by radio halos, from cluster core (solid 
lines) to 1 Mpc distance from the cluster center (dashed 
lines); (Vl/c b ) 2 =0.22 is assumed in both cases. Figure 4 
highlights the effect of decreasing radiative and Coulomb 
losses on the spectrum of the reaccelerated electrons : the 
boosting of the particle spectrum increases in the exter- 
nal (Mpc-distance) regions because Coulomb losses are less 
severe (which allows the reacceleration of the bulk of sec- 
ondary electrons) and because synchrotron losses become 
sub-dominant with respect to inverse Compton losses due 
to the scattering of the CMB photons. 



r kk 



4.3 The non-thermal spectrum of turbulent 
galaxy clusters 

In this Section we calculate the non thermal spectrum from 
galaxy clusters assuming reacceleration of primary protons 
and of their secondary products by compressible MHD tur- 
bulence in the IGM. 
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Figure 5. The radio (left) and high energy (right) emitted spectra from the Coma cluster (note that the SZ decrement at higher radio 
frequencies is not included, see e.g. Donnert et al. 2010b). Upper panels show our reference model (see text), lower panels refer to the 
case of a cluster magnetic field 2.5 times smaller than that in our reference model (for a better comparison red dot-dashed lines in the 
right panel show the IC and gamma-ray emission expected in our reference model). Solid lines show the emitted spectra obtained after 6.5 
xlO 8 yrs of reacceleration (upper and lower lines in the synchrotron and gamma-ray spectra refer to the emission within 4.5 and 3 core 
radii, respectively, while X— rays are calculated within 4.5 core radii), dashed lines show the emitted spectra obtained by assuming that 
MHD turbulence is dissipated in the cluster. Dotted lines show the emitted spectra at intermediate stages of reacceleration : after 2 and 
4xl0 8 yrs, from bottom to top respectively. Radio data-points are taken from Thierbach et al (2003), the hard X-ray upper limit is from 
Suzaku observations (Wik et al 2008) and the EGRET and FERMI upper— limits are taken from Reimer et al (2003) and Ackermann et al 
(2010), respectively. Calculations assume (V L /csf ~ 0.18 for our reference model and 0.2 for the lower panels (L = 300 kpc is assumed 
in all calculations). For reference, the energy of the population of relativistic protons "measured" after 0.65 Gyr of reacceleration (solid 
lines) is 0.035 and 0.25 times that of the thermal cluster in our reference model and lower panels, respectively. 
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Figure 4. The time evolution of the spectrum of electrons in 
two different environments: in the clusters central regions (solid 
lines), assuming n th = 2.64 X 10~ 3 cm — 3 and B = 4.69^tG, 
and in the clusters external regions (red dashed lines), assum- 
ing n th = 0.56 X 10 _3 cm -3 and B = 2.16/iG. From bottom 
to top, calculations are shown after 0.3, 2, 4 and 6.5 xlO 8 yrs 
of reacceleration. Calculations are obtained at z = 0, assuming 
(Vl/cs) 2 = 0.22, with injection scale L a = 300 kpc, T = 10 8 K. 



4-3.1 A toy model for the Coma cluster 

The Coma cluster hosts the best studied, prototype, gi- 
ant radio halo (Willson 1970; Giovannini et al 1993). Con- 
straints on the high energy emission from the Coma cluster 
are presently available from hard X-ray (BeppoSAX, Fusco- 
Femiano et al 1999, 2004, Rossetti & Molendi 2004; RXTE, 
Rephaeli et al 1999; INTEGRAL, Eckert et al 2007, Lu- 
tovinov et al 2008; Suzaku, Wik et al 2009; Swift-BAT, 
Ajello et al 2009) and gamma ray (EGRET, Reimer et al 
2003; FERMI, Ackermann et al 2010; HESS, Aharonian et 
al 2009b; VERITAS, Perkins 2008) observations. It is thus 
a natural first step to compare our model expectations with 
the spectral energy distribution (SED) of the non thermal 
emission from the Coma cluster. 

We assume a simplified model for the thermal gas distri- 
bution in the Coma cluster that is anchored to the observed 
beta-model profile (e.g. Briel et al 1992, considering ACDM 
cosmology). The gas temperature, ksT ~ 8 keV (David et 
al 1993), is assumed to be constant on Mpc-scale. 
The magnetic field in the Coma cluster, and its spatial dis- 
tribution, is a crucial ingredient in our modeling of the syn- 
chrotron (radio halo) properties. We adopt the recent results 
by Bonafede et al (2010) that carried out a detailed analy- 
sis of the Rotation Measures of a sample of extended cluster 
radio galaxies in the Coma cluster. The best fit to their data 
implies a magnetic field B(r) oc sjrith with a central value 
B(0) ~ 5/iG. 

Once the thermal gas and magnetic field properties are 
anchored to present data, the free parameters in our calcu- 
lations are (i) the energy density (and spatial distribution) 



of relativistic protons (at time = 0)++, and (ii) the injection 
rate (and scale) of compressible turbulence in the IGM, I . 
In our reference toy model we assume (i) a flat spatial distri- 
bution of relativistic protons on the halo-scale, r ~ 3 r c , and 
a constant ratio between relativistic and thermal particles 
energy densities at larger distances, and (ii) a specific in- 
jection rate of turbulence, I /p, constant on the halo-scale. 
Assumption (i) is motivated by the fact that, similarly to 
other giant halos, the Coma radio halo has a very broad 
synchrotron-brightness distribution (Govoni et al 2001; see 
also discussion in Cassano et al 2007 and Donnert et al 
2010a) implying a very broad spatial distribution of rela- 
tivistic protons on the halo-scale. Assumption (ii) is mo- 
tivated by hydrodynamical and MHD cosmological simula- 
tions that found very extended turbulent regions in simu- 
lated galaxy clusters (Sunyaev et al 2003; Dolag et al. 2005; 
Vazza et al 2009a; Paul et al 2010). 

In Figure 5 (upper panels) we show the expected SED 
emitted by the Coma cluster (from r < 3 — 4.5 r c ) assuming 
our reference model (see caption). We find that the radio 
spectrum of the Coma radio halo can be well reproduced by 
assuming a total energy content of relativistic protons, on 
the halo-scale, Ecu ~ 3.5% of the thermal gas (see caption 
for details). The presence of a break in the spectrum of the 
Coma radio halo at higher frequencies has been interpreted 
as a signature of turbulent acceleration^ , since it implies a 
corresponding break in the spectrum of the emitting elec- 
trons at the energy where turbulent acceleration is balanced 
by radiative losses (e.g. Schlickeiser et al 1987; Brunetti et 
al 2001; Petrosian 2001). Figure 5 demonstrates that also 
models considering the reacceleration of secondary particles 
can explain the steepening of the Coma radio halo, although 
a tail at higher frequencies due to freshly-injected secondary 
electrons shows up in the (emitted) spectrum^ . For a con- 
sistency check, in Figure 6 we also show that, although very 
simplified, our reference model allows us to (roughly) repro- 
duce the observed radio brightness profile of the Coma halo 
(although the predicted profile is still slightly steeper than 
the observed one). Such expected-broad synchrotron profile 
is due to the combination of the flat spatial distribution of 
relativistic protons with the increasing efficiency of reaccel- 
eration of secondary electrons at larger distances from the 
cluster centre (see Fig. 4). 

The gamma ray emission from the cluster (from r < 
3 — 4.5 r c , see caption) is dominated by the decay process of 
7T° (Fig. 5) and is expected at about 10 percent level of the 
present gamma-ray upper limits (FERMI, Ackermann et al 
2010). In our model the energy content of relativistic protons 
that is necessary to reproduce the observed luminosity and 
brightness profile of the Coma radio halo is much smaller 
than that from calculations based on pure hadronic models 
that assume a similar magnetic field in the Coma cluster 
(Donnert et al 2010b), implying also a much smaller gamma- 
ray luminosity. 



M we assume an initial spectrum of protons N p (p) = K v p 2,6 

more recently Donnert et al. (2010b) have shown that the 
observed steepening cannot be due to the SZ decrement due to 
the hot gas in the central Mpc region of the cluster 

a similar conclusion comes from calculations of reacceleration 
of secondary particles in the Alfvenic case (Brunetti & Blasi 2005) 
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Figure 6. Radio brightness profile of the Coma radio halo (in ar- 
bitrary units). The observed profile (tick-blue lines and crosses) is 
taken from Govoni et al (2001) observations. Solid lines show the 
brightness profile at 330 MHz predicted by our reference model 
(red line) and by the model with lower magnetic field (green line). 
Dashed lines refer to the brightness profiles of both models (same 
color code) considering the case where turbulence is dissipated. 



For seek of completeness, in Figure 5 (lower panels) we show 
the expected SED of the Coma cluster by assuming a mag- 
netic field in the cluster 2.5 times smaller than that inferred 
from Bonafede et al.(2010) (see caption). In this case the 
properties of the Coma radio halo can be reproduced by as- 
suming an energy content of relativistic protons Ecu ~ 25% 
of the thermal energy, i.e. 7-8 times larger than that of our 
reference model, and the expected gamma ray emission in- 
creases significantly. Interestingly, in this case there would 
be a chance to detect the Coma cluster in the next years 
with the FERMI telescope. 

Remarkably, on the other way round Figure 5 (lower panels) 
demonstrates the importance of present gamma ray upper 
limits: in the turbulent reacceleration picture, the magnetic 
field in the central regions of the Coma cluster cannot be sig- 
nificantly smaller than about 2 pG, provided that hadronic 
collisions are the main source of the (seed) electrons in the 
cluster volume. 

In the case of our reference model, the inverse Compton 
emission in the hard X-ray band is expected at a ~percent 
level than the present Suzaku upper limit and a detection 
of the Coma cluster will be challenging even with the fu- 
ture hard X-ray experiments (e.g. ASTRO-H, Nustar). On 
the other hand, by assuming a magnetic field in the clus- 
ter 2.5 times smaller than that inferred from Bonafede et 
al.(2010) the expected inverse Compton emission increases 
significantly and there would be a chance to detect the Coma 
cluster with future hard X-ray telescopes. However, by as- 
suming a magnetic field substantially smaller than that of 
our reference model, we find that the predicted radio profile 
is much steeper than the observed one (Fig. 6). 



4-3.2 Transient and long-living spectral components in 
galaxy clusters 

In the context of our model the non-thermal emission from 
galaxy clusters is a mixture of two main spectral components 
: a long-living one that is emitted by the chain of secondary 
particles continuously generated by collisions between ther- 
mal and long-living (several Gyrs) relativistic protons, and 
a transient amplification of the SED that appears when rel- 
ativistic particles are reaccelerated by the MHD turbulence 
generated (and then dissipated) in connection with cluster 
mergers (see also Brunetti et al 2009a); in this scenario the 
idea is that the transient synchrotron component generates 
the observed radio halos. 

The non-thermal spectrum of dynamically relaxed galaxy 
clusters should be mainly due to the long-living spectral 
component. An example of this spectral component is shown 
in Figure 5 (dashed lines) where we report the synchrotron, 
inverse Compton and n° emission calculated by assuming 
that MHD turbulence in the Coma cluster is dissipated. 
The main effect of turbulent-reacceleration is to produce a 
amplification of the synchrotron (radio), inverse Compton 
(hard X) and 7r° (gamma) emission in merging clusters (Fig. 
5), while the effect of the dissipation of MHD turbulence in 
those clusters that become more relaxed is to suppress (only) 

the radio and the hard X-ray emissiorJUL Consequently in 
our model we expect a tight correlation between clusters 
dynamics and the radio (radio halos) and hard X-ray prop- 
erties of galaxy clusters, but do not expect a tight correlation 
between gamma rays and cluster dynamics. 

An important issue is the evolution of the radio halos in 
connection with cluster mergers and the difference between 
the radio properties of merging and relaxed clusters. 
In the "classical" scenario, where turbulence reaccelerates 
primary electrons, the boosting of the synchrotron and hard 
X-ray emission during cluster mergers can be extremely 
large (eg. Fig. 5 in Brunetti et al 2009b; Fig. 6 in Cas- 
sano 2010). Contrary to that, in our model the presence of 
relativistic protons in the IGM and the chain of secondary- 
particles decay induced by pp collisions put stringent con- 
straints to the radio-evolution of galaxy clusters, since 
long-living spectral components are generated for about a 
Hubble time. Our calculations exploit this point and un- 
veil the interplay between transient and long-living spectral 
components. In our model the boosting of the radio and 
hard X-ray luminosities due to reacceleration is fairly con- 
strained, since it is essentially due to the reacceleration of 
secondary electrons from the energies where the injection 
spectrum of secondary electrons peaks (or from the energy 
where reacceleration is stronger than Coulomb losses) to 
the energies necessary to electrons to radiate synchrotron 
emission in the radio band (provided reacceleration is suf- 
ficient to accelerate particles up to these energies). This 
is SP/P < (EQ ± (E)) El /(EQ ± (E)) Eh , where E t and E h 
are about 100 MeV and 1 GeV respectively, that implies 
5P/P ~ 10 — 15 for typical values of the slope of the proton 
spectrum, s = 2.4 — 2.8. 

The predicted level of amplification and suppression of radio 



HI Also the radial synchrotron profile in non-turbulent clusters 
is steeper than that in turbulent clusters (Fig. 6) 
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emission in galaxy clusters and its connection with cluster 
dynamics can be constrained by radio and X-ray observa- 
tions. The recent radio follow up of a complete X-ray sample 
of galaxy clusters, the " GMRT Radio Halo Survey" ( Venturi 
et al. 2007, 2008), leads to the discovery of a bi-modal be- 
haviour of the clusters radio properties (Brunetti et al 2007, 
2009b): radio halos are found in the most disturbed systems 
and trace the Pi. 4 — Lx correlation, on the other hand the 
majority of clusters are less disturbed and "radio quiet", 
with the limits to their radio luminosities about 10 times 
smaller than the radio luminosities of halos. As a matter of 
fact, our reacceleration model predicts that the level of the 
long- living synchrotron radiation in "radio quiet" clusters 
is similar to that of the upper limits from current GMRT 
radio observations, on the other hand the "classical" reac- 
celeration scenario, where secondaries do not play a role, 
would predict a larger gap between radio halos and "radio 
quiet" clusters. Consequently future (deeper) radio observa- 
tions may test the role of relativistic protons and secondary 
particles. 

A second important issue is the possibility to detect 
gamma ray emission from galaxy clusters. In our reacceler- 
ation model gamma ray emission is a spectral component 
common in galaxy clusters, it essentially depends on the 
integrated merging history of clusters rather than to the dy- 
namical status of a given cluster at the epoch of the observa- 
tion; this is similar to classical hadronic models (eg., Blasi et 
al 2007 for review) and indeed our model evolves into these 
models when MHD turbulence in the IGM is dissipated. 
In our "reference" model-configuration we expect that 
gamma ray emission from galaxy clusters is difficult to de- 
tect in next years with present facilities. However, assuming 
a smaller magnetic field in the IGM, a larger number of 
relativistic protons must be assumed to explain radio halos 
and there is a chance to detect gamma rays from clusters 
in next years. As in the case of hadronic models, we expect 
that gamma ray emission scales with thermal properties as 

L 7 oc p 2 IGM MT{^) (30) 

where (...) indicates emission weighted quantities in the clus- 
ter volume. In the most simple situation where we assume 
that (^ a ) does not depend on cluster mass, Eq.30 implies a 
quasi linear scaling between gamma and (thermal) X-ray lu- 
minosities of galaxy clusters, L 1 oc L x 3 * * * (see also Pfrom- 
mer 2008 and Donnert et al 2010a for more detailed calcu- 
lations based on cosmological simulations). 



5 DISCUSSION AND CONCLUSIONS 

If MHD turbulence is generated in galaxy clusters during 
massive merger events, relativistic protons can be reacceler- 
ated enhancing the rate of generation of secondary particles 
in the IGM. Consequently, for a self-consistent treatment of 
the process, protons and their secondary products must be 
taken into account in the calculations of particle reaccelera- 
tion by MHD turbulence in galaxy clusters. 

* * * Here we have considered Lx oc p^ GM M\/T, and assumed 
the empirical scaling Lx oc T 2 - 98 (eg. Reiprich 2001). 



In this paper we have assumed the picture of MHD turbu- 
lence in galaxy clusters as described in Brunetti & Lazarian 
(2007) to calculate the acceleration of primary protons and 
secondary particles by compressible MHD turbulence. 

Upper limits to the gamma ray emission from galaxy 
clusters and radio upper limits to the cluster-scale synchro- 
ton emission in clusters without radio halos constrain the 
energy density of relativistic protons in the central Mpc re- 
gions of galaxy clusters to less than a few percent of that 
of the IGM (Reimer et al 2004; Brunetti et al 2007, 2008; 
Aharonian et al 2009a,b; Ackermann et al 2010). 

We have shown that a population of relativistic protons 
consistent with the above limits may be sufficient to gener- 
ate the observed radio halos in merging clusters (including 
their observed brightness profiles and spectra) via the reac- 
celeration of their secondary electrons by compressible MHD 
turbulence, provided that the spatial distribution of the rel- 
ativistic protons is relatively flat and that the magnetic field 
in galaxy clusters is at the level derived from Rotation Mea- 
surements. 

The main consequence of this theoretical scenario is that the 
non-thermal SED of galaxy clusters is given by the interplay 
of a transient component, generated by the reacceleration of 
secondary particles by MHD turbulence during cluster merg- 
ers (that generates radio halos), and a long-living compo- 
nent, that is generated by the secondary particles that are 
continuously injected by pp collisions in the IGM. In this 
case we expect a tight connection between radio halos and 
cluster mergers, as well as an amplification of the level of 
hard X-ray emission in merging clusters, while we expect no 
tight correlation between gamma rays and cluster dynamics. 
At the same time, we also expect that diffuse radio emission, 
due to secondary particles, must be common in galaxy clus- 
ters at a level that is about one order of magnitude below 
that of nowadays observed radio halos, i.e. at the same level 
of the upper limits derived by present radio observations of 
"radio quiet" clusters. This expectation can be tested by 
future (deeper) radio observations of "radio quiet" clusters. 
The level of gamma ray emission from nearby, massive, 
galaxy clusters is expected at about 10% of the level of 
present upper limits, assuming a magnetic field strength in 
these clusters in line with present studies of RM. This im- 
plies that only future telescopes (e.g. CTA) may lead to the 
detection of galaxy clusters in the gamma ray band. On the 
other hand, if the magnetic field in galaxy clusters is smaller, 
detection of galaxy clusters with FERMI could be possible 
in next years. In this respect, in the context of our model, 
present FERMI upper limits for the Coma cluster already 
provide a limit to the central value of the magnetic field 
B(0) > 1 - 2^G. 

5.1 Model simplifications and future steps 

In our paper we focus on the role played by compress- 
ible MHD turbulence, essentially the fast modes. Our ed- 
ucated guess, motivated in Sect. 2 and in Brunetti & Lazar- 
ian (2007), is that these modes are the most relevant for 
the reacceleration of relativistic particles in galaxy clusters. 
These modes were also identified as major scattering agent 
for Galactic cosmic rays (Yan & Lazarian 2002, 2004). 
The Alfven modes and slow modes are not efficient for scat- 
tering if the turbulent energy is being injected at large 
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scales (Chandran 2000, Yan & Lazarian 2002). Their in- 
efficiency stems from the both the spectra being steep in 
terms of parallel perturbations as well as fluctuations be- 
ing very anisotropic (Goldreich & Sridhar 1995, Lazarian 
& Vishniac 1999, Cho & Vishniac 2000, Maron & Goldre- 
ich 2001, Cho, Lazarian & Vishniac 2002). If Alfven modes 
are injected by instabilities they may have radically differ- 
ent properties from the modes of the large scale cascade. 
For instance, Alfven modes arising from particle streaming 
have slab structure. Such waves efficiently interact with en- 
ergetic particles. Other instabilities, e.g. gyroresonance one, 
can produce slab waves (see Gary 1993). The instabilities 
producing slab Alfven modes may be induced by large scale 
compressible turbulence (see Lazarian & Beresnyak 2006). 
Indeed it is worth mentioning that the mode composition at 
smaller scales, I << I a, in the IGM could becomes rather 
complex (e.g. Kato 1968; Eilek & Henriksen 1984). 

To what extend our calculations are accurate depends 
on our understanding of properties of IGM turbulence. We 
assumed that the damping of fast modes with thermal par- 
ticles can be obtained considering a collisionless plasma. 
At the same time, it can be argued (see e.g., Lazarian et 
al. 2010) that the degree of collisionality of astrophysical 
plasmas can be underestimated if only Coloumb collisions 
are taken into account, as particles in plasmas can inter- 
act through the mediation of the perturbed magnetic fields. 
The inevitable conclusion is that the collisionless formulae 
describing damping of these modes should be only applied 
to scales less than the mean free path, which is much shorter 
than the Coulomb mean free path. As a result, fast modes 
should be substantially less damped and be present on the 
scales which are much shorter than the earlier estimates, in- 
cluding those in this paper. The consequence of this is that a 
more appreciable portion of energy gets available for the ac- 
celeration of cosmic rays. Therefore, our present calculations 
may underestimate the efficiency of cosmic ray acceleration 
by turbulence. We plan to address the self-consistent prob- 
lem elsewhere. 

Following Brunetti & Lazarian (2007) in our paper we 
adopt quasi-linear-theory (henceforth QLT) to calculate 
particle acceleration by fast modes. In Brunetti & Lazarian 
(2007) we indeed have shown that for fast modes in the IGM 
it is (cj) >> (r), where (...) indicates angle-averaged quan- 
tities, that provides some justification to the use of QLT. 
More recent studies in Yan, Lazarian & Petrosian (2008) pre- 
sented an approach to particle acceleration that generalizes 
the QLT and allows to take into account the effect of large 
scale variations of magnetic field. As numerical simulations 
which use the data from the actual MHD turbulence simu- 
lations (Beresnyak, Yan, Lazarian 2010) support the theory, 
we believe that, in future, going beyond the standard QLT 
approach may provide a more accurate description of the 
particle acceleration process, although we do not expect a 
large change (see Yan et al 2008). 

Finally a unavoidable simplification in our semi- 
analytical calculations is that turbulence is homogeneous, 
in space and time, on the radio halo volume, during cluster 
mergers. On the other hand, cosmological numerical sim- 
ulations of galaxy clusters show a more complex situation 
where intermittent and patchy large-scale turbulent motions 
are generated during multiple collisions between galaxy clus- 
ters (e.g. Vazza et al 2009a; Paul et al 2010). We believe 



that to be more realistic one may need to vary the intensity 
of turbulence driving to describe different parts of galaxy 
clusters. Our idealized calculations of particle acceleration 
and evolution of compressible turbulence in the IGM pro- 
vides a first step, implementing our formalism in detailed 
cluster-simulations is necessary to obtain a more reliable 
description of the morphology and spectral distribution of 
the non-thermal emission and of the connection between 
radio halos and cluster mergers. 

5.2 A comparison with the Alfvenic approach 

Earlier calculations of turbulent acceleration of primary pro- 
tons and secondary electrons in galaxy clusters focus on the 
reacceleration by Alfven modes (Brunetti & Blasi 2005). 
These studies first provided a description of the expected 
transient and long-living spectral components in galaxy 
clusters (e.g., Brunetti et al 2009a). 

It is well known that the damping of these modes is mainly 
due to the interaction (gyro-resonance) with relativistic par- 
ticles, that provides the main motivation to explore this 
possibility. On the other hand, as already mentioned, a 
complication of this approach is the anisotropy of Alfven 
modes that develops when turbulence cascades from larger 
to smaller scales. Consequently in Alfvenic models the injec- 
tion of Alfven modes "directly" at small scales, i.e. compa- 
rable with the gyroradius of high energy particles, must be 
postulated, in which case it is also difficult to derive a over- 
all picture connecting clusters mergers and the generation 
of these modes at such small scales. 

A second issue in the modeling of Alfvenic acceleration of 
relativistic protons and of their secondary products is that 
the primary and secondary particles interact with Alfven 
modes with different scales. This is because secondary elec- 
trons of energy E e are mostly generated by collisions be- 
tween higher energy protons, E p ~ 10-50 xE e , and ther- 
mal targets that implies that these secondaries interact with 
modes on scales 10-50 times smaller than those of the par- 
ent primary protons (by consider the gyroresonant condi- 
tions k ~ eB/(cp cos#)). Consequently the ratio between 
the transient and long-living components in these models 
depends also on the spectrum of Alfven modes. This implies 
a larger degree of freedom in Alfvenic models with respect 
to the more straightforward case treated in our paper where 
compressible MHD turbulence interact with relativistic par- 
ticles. 



5.3 Sources of primary electrons in a turbulent 
IGM 

In general, if relativistic protons and electrons are present 
in the IGM, we expect that the MHD turbulence, generated 
during cluster mergers, would reaccelerate both these par- 
ticles. If we assume an efficient confinement of cosmic rays 
in galaxy clusters, the unavoidable consequence of this sce- 
nario is that also the energy density of secondary products, 
due to pp collisions, should increase and consequently their 
contribution to the non-thermal clusters spectrum. 
In our paper we have addressed this problem under the most 
extreme (and simplified) condition where only protons and 
their secondaries are present in the IGM; we note that un- 
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der these conditions the ratio of the energy densities of rel- 
ativistic protons and of the emitting electrons in the IGM 
is maximized. In this context, we have shown that, assum- 
ing turbulent reacceleration at the level necessary to explain 
radio halos, the radio emission generated by secondary parti- 
cles when turbulence is dissipated is consistent with present 
upper limits to the diffuse radio luminosity of clusters with- 
out radio halos. Also we have shown that the gamma-ray 
emission from the decay of secondary neultral-pions, that 
should be common in galaxy clusters and not tightly con- 
nected with their dynamical status, is expected at a level 
much smaller than that constrained from present upper lim- 
its with FERMI. 

At the same time, however, under these conditions our cal- 
culations suggest that the difference between the cluster- 
scale radio emission of "turbulent" (merging) and "non- 
turbulent" (relaxed) clusters cannot be larger than about a 
factor ten, and that we expect to detect diffuse radio emis- 
sion potentially in "all" massive clusters as soon as much 
deeper observations of clusters that are presently defined 
"radio-quiet" will become available. 

On the other hand primary relativistic electrons should 
be present, at some level, in the IGM. It is well known 
that active radio galaxies may fill large volumes in the IGM 
with relativistic plasma, relativistic electrons age rapidly 
but they can be accumulated at hundred MeV energies for 
longer times, especially in the external regions (Fig.l). Other 
sources of relativistic primary electrons that are usually con- 
sidered in the literature are Galactic Winds and Starburst 
galaxies (e.g., Volk & Atoyan 1999) and shock waves (e.g., 
Ensslin et al 1998; Sarazin 1999; Ryu et al 2003; Pfrommer et 
al. 2006; Skillman et al 2008; Vazza et al 2009b), evidence 
for the latter process come from the observations of radio 
relics that indeed suggest a connection between shocks and 
electron acceleration (or reacceleration) in the IGM (e.g., 
Markevitch et al 2005 and Giacintucci et al 2008 for obser- 
vations of a shocks-relics connection in clusters). 
In addition to these processes, another possibility that re- 
quires more attention is that the primary electrons can be 
"created" in the IGM. Magnetic reconnection presents a 
natural way for doing this, as electrons bounces back and 
forth between converging magnetic fluxes can gain energy 
through the first order Fermi acceleration (Gouveia dal Pino 
& Lazarian 2003). Although it is generally believed that re- 
connection is a slow process, potentially turbulence can sig- 
nificantly enhance the reconnection rate (Lazarian & Vish- 
niak 1999). Exploring this mechanism and its interplay with 
particle reacceleration by MHD turbulence in the IGM may 
open new perspective in our understanding of non-thermal 
cluster-scale emission and of its connection with cluster 
mergers. We aim to discuss in detail this point in a future 
paper, yet for completeness here we dedicate an Appendix 
to turbulent reconnection. 

The presence of primary electrons in the IGM may sig- 
nificantly affect the SED of galaxy clusters and its evolution. 
In particular if the number density of primary electrons is 
comparable to (or larger than) that of secondary particles, a 
smaller number of relativistic protons is required by the our 
model to match the observed spectrum of radio halos. This 
has an impact on the expected gamma ray emission from 
galaxy clusters that would be consequently smaller than that 
expected from the calculations presented in Sect. 4. 3. A sec- 



ond important point is that if a population of reaccelerated 
primary electrons significantly contributes to the observed 
radio halo emission, the expected ratio between transient 
and long living components in the radio and hard X-ray 
bands increases with respect to that derived in Sect. 4. 3. 
As a matter of fact, assuming that primary electrons are 
dominant with respect to secondaries our picture evolves 
into the "classical" reacceleration model (eg., Brunetti et al 
2001, 2004; Petrosian 2001; Cassano & Brunetti 2005) where 
indeed the luminosity of giant radio halos is suppressed by 
several orders of magnitude when galaxy clusters evolve into 
relaxed systems (e.g. Brunetti et al 2009b; Cassano 2010). 
Consequently, we believe that future gamma ray and radio 
observations will be crucial also to constrain the ratio of 
primary and secondary electrons in the IGM. 
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7 APPENDIX: MAGNETIC FIELD 
RECONNECTION AS SOURCE OF 
PRIMARY ELECTRONS IN THE IGM 

Magnetic reconnection is an ubiquitous process in magne- 
tized flows, expected to happen when magnetic fields of 
non-parallel direction get into contact. However, the text- 
book processes of magnetic reconnection seem to fall short 
of providing the desired solution. The Sweet-Parker recon- 
nection (Sweet 1958, Parker 1957) is too slow. The speed 
of the reconnection scales as Alfven velocity va times the 
inverse value of the square root of the Lunquist number 
S* 1 / 2 = (Lva/i]) 1 ^ 2 , where L is the size of the magnetic re- 
gions and r) is magnetic diffusivity. As Lunquist number can 
be for the IGM 10 15 and higher, it is clear that the Sweet- 
Parker reconnection can handle only a negligible fraction 
of magnetic flux in the Hubble time. Petscheck proposed a 
solution where magnetic fields get into contact at a sharp 
angle (Petscheck 1964). 

In the Petscheck model no efficient acceleration of par- 
ticles is expected at reconnection sites. Indeed the tradi- 
tional processes of acceleration which rely on the electric 
field in the reconnection region are inefficient, as the recon- 
nection region is too small and releases only a small portion 
of magnetic energy. Slow shocks predicted in the Petscheck 
model are likely to be inefficient for the particle acceleration 
(see discussion in Beresnyak, Jones & Lazarian 2009). In 
addition, observational data does not support the X-point 
reconnection predicted in the Petscheck model either (see 
Ciaravella k, Raymond 2008). 

A shortcoming of many discussions of magnetic recon- 
nection is that the traditional setup does not include ubiq- 
uitous pre-existing astrophysical turbulence. As turbulence 
radically changes many astrophysical processes, the influ- 
ence of turbulence on reconnection has attracted the atten- 
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Figure 7. Cosmic rays entrained over a shrinking loop of recon- 
nected magnetic flux bounce between points A and B and gain 
energy (from Lazarian 2005). 



tion of researchers for a long time (see Speizer 1970; Math- 
aeus & Lamkin 1985, 1986; Strauss 1988). 
A new approach to the effects of turbulence was adopted in 
Lazarian & Vishniac (1999, henceforth LV99). This model 
predicts reconnection speeds close to the turbulent velocity 
in the fluid. More precisely, assuming isotropically driven 
turbulence characterized by an injection scale, I, smaller 
than the current sheet length L, LV99 obtained : 

V rec ^v A (l/L) 1/2 (V,/v A ) 2 , (31) 

where the turbulent injection velocity Vi is assumed to 
be less than va- If L < I, the first factor in Eq. (31) should 
be changed to (L/l) 1 ^ 2 (LV99). If turbulent injection veloc- 
ity is larger than va the reconnection happens at the Alfven 
speed for L > Ia- For L < Ia a factor (L/Ia) 1 '" 2 should 
substitute the factor of (l/L) 1 '" 2 in Eq. (31). Physically this 
reflects the fact that at sufficiently small scales magnetic 
field energy dominate the kinetic energy and the magnetic 
field lines get only weakly perturbed by turbulent motions. 
Figure 7 provides the simplest realization of the acceleration 
of particles within the reconnection region expected within 
LV99 model. As a particle bounces back and forth between 
converging magnetic fluxes, it gains energy through the first 
order Fermi acceleration (de Gouveia dal Pino & Lazarian 
2003, 2005; see also Lazarian 2005). The first order accelera- 
tion of particles entrained on contracting magnetic loop can 
be understood from the Liouville theorem, i.e the preser- 
vation of the phase volume which includes the spatial and 
momentum coordinates. As in the process of reconnection 
the magnetic tubes are contracting and the configuration 
space presented by magnetic field shrinks, the regular in- 
crease of the particle's energies is expected. The requirement 
for the process to proceed efficiently is to keep the acceler- 
ated particles within the contracting magnetic loop. This 
introduces limitations on the particle diffusivities perpen- 
dicular to magnetic field direction. Thus high perpendicular 
diffusion of particles may decouple them from the magnetic 
field. Indeed, it is easy to see that while the particles within 
a magnetic flux rope depicted in Figure 7 bounce back and 
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forth between the converging mirrors and get accelerated, if 
these particles leave the flux rope fast, they may start bounc- 
ing between the magnetic fields of different flux ropes which 
may sometimes decrease their energy. Thus it is important 
that the particle diffusion parallel and perpendicular mag- 
netic field stays different. Particle anisotropy which arises 
from particle preferentially getting acceleration in terms of 
the parallel momentum may also be important. The energy 
spectrum was derived in GL03 : 

N(E)dE = C ■ E~ 5/2 dE, (32) 

In Drake et al. (2006) this idea was enriched by taking 
into account the backreaction of particles. 
We believe that such an acceleration can be present in 1GM. 
Two dimensional numerical simulations of acceleration can 
be found in Drake et al. (2006, 20f0) and first three dimen- 
sional simulations were presented in Lazarian et al. (20f0). 
In typical IGM the small scale reconnection of turbulent field 
should be collisionless and during this collisionless with elec- 
trons and ions decoupled. This should enable the accelera- 
tion of electrons. The idea of particle acceleration in recon- 
nection regions was recently applied to describe anomalous 
cosmic ray acceleration in Heliosphere (Lazarian &: Opher 
2009, Drake et al.2010) and acceleration of cosmic rays in 
heliotail (Lazarian & Desiati 2010). In IGM magnetic tur- 
bulence creates magnetic reversals and therefore we expect 
to see additional sources of energetic electrons. In terms of 
the model that we discuss the electron acceleration in mag- 
netic reconnection sites should increase the energy density of 
electrons and modify their spectrum. We shall discuss these 
effects elsewhere. 
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